Gaussian-state theory of two-photon imaging 
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Biphoton states of signal and idler fields — obtained from spontaneous parametric downconversion 
(SPDC) in tlie low-briglitness, low-flux regime — have been utilized in several quantum imaging 
configurations to exceed the resolution performance of conventional imagers that employ coherent- 
state or thermal light. Recent work — using the full Gaussian-state description of SPDC — has shown 
that the same resolution performance seen in quantum optical coherence tomography and the same 
imaging characteristics found in quantum ghost imaging can be realized by classical-state imagers 
that make use of phase-sensitive cross correlations. This paper extends the Gaussian-state analysis 
to two additional biphoton-state quantum imaging scenarios; far field diffraction-pattern imaging; 
and broadband thin-lens imaging. It is shown that the spatial resolution behavior in both cases is 
controlled by the nonzero phase-sensitive cross correlation between the signal and idler fields. Thus, 
the same resolution can be achieved in these two configurations with classical-state signal and idler 
fields possessing a nonzero phase-sensitive cross correlation. 

PACS numbers: 42.30.Va, 42.50.Ar, 42.50.Dv 



I. INTRODUCTION 

Spontaneous parametric downconversion (SPDC), 
with a continuous-wave nondepleting pump, produces 
signal and idler fields that are in a maximally-entangled, 
zero- mean jointly Gaussian state [J, 01 ■ When the 
downconverter is operated in its low-brightness, low- 
flux regime, so that at most one signal-idler photon 
pair is emitted during an observation interval, this state 
reduces to the superposition of a predominant multi- 
mode vacuum state plus a weak biphoton component Q . 
Biphoton illumination is at the heart of several quantum 
imaging configurations, including quantum optical coher- 
ence tomography (OCT) [1], quantum ghost imaging 5|, 
quantum holography [6| and quantum lithography 7|. 
These systems offer performance advantages over con- 
ventional optical imagers, which employ coherent-state 
or thermal-state sources, that have traditionally been as- 
cribed to the entanglement between the biphoton's signal 
and idler components. We have shown — using Gaussian- 
state analysis and phase-sensitive coherence theory — that 
the advantages in quantum OCT and ghost imaging pre- 
dominantly stem from the phase-sensitive cross correla- 
tion between the signal and idler fields, rather than their 
entanglement per se [j, Furthermore, because a pair 
of classical-state fields TO'] can also have nonzero phase- 
sensitive cross correlation, most of the advantages seen 
in these biphoton-state imagers are also attainable with 
classical phase-sensitive sources, but conceivably (and 
conveniently) at much higher photon flux and without 
the need for single-photon counters. 

In this paper we will study Fourier-plane and thin-lens 
imaging of a transmission mask using jointly-Gaussian 
source states. These states encompass the biphoton 
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state, thermal states (used in conventional low-coherence 
imaging) and coherent states (used in conventional co- 
herent imaging). Previous theoretical and experimen- 
tal work has shown that biphoton-state illumination of a 
transmission mask yields a far-fleld diffraction (Fourier- 
plane) pattern that is factor-of-two compressed rela- 
tive to that produced by coherent plane-wave illumina- 
tion of the mask In addition, it has been claimed 
that imaging a transmission mask using broadband, 
spatially-incoherent biphoton-state illumination and a 
finite-diameter thin lens will yield a point-spread func- 
tion that is a factor of two narrower than that obtained 
with quasimonochromatic spatially-incoherent thermal- 
state illumination fTT| . As noted above for quantum OCT 
and ghost imaging, these biphoton-state resolution en- 
hancements in Fourier-plane and thin-lens imaging have 
been ascribed to signal-idler entanglement. We will use 
Gaussian-state analysis, together with phase-sensitive co- 
herence theory, to develop a unified — and generalized — 
understanding of the classical and quantum regimes of 
these imaging configurations. In particular, we will show 
that, once again, the benefits ascribed to entanglement 
are in fact due to phase-sensitive coherence at the source, 
and therefore are obtainable with classical-state sources 
with phase-sensitive coherence. Furthermore, we will 
show that the narrowing of the point-spread function in 
the thin-lens imaging configuration is much less than a 
factor of two even at the theoretical upper limit of SPDC 
bandwidth, and that it is marginal for typical SPDC 
bandwidths. 

The remainder of this paper is organized as follows. 
In Sec. [ni we present the phase-insensitive and phase- 
sensitive coherence theory results that are relevant to 
Fourier-plane imaging. In Sec. IIIII we use the Sec. |TT] 
results for quasimonchromatic light to determine the far- 
field diffraction properties of Gaussian-state source fields 
after they have illuminated a source-plane transmission 
mask. Section HVl presents our analysis of thin- lens imag- 
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ing, in which we determine the effect of the source's 
bandwidth on the point-spread functions resulting from 
phase-insensitive and phase-sensitive Gaussian-state il- 
lumination. In doing so we pay particular attention to 
frequency-dependent propagation effects that did not en- 
ter into our quasimonochromatic treatment of Fourier- 
plane imaging. Finally, in Sec. |Vl we summarize the 
results we have obtained for the Fourier-plane and thin- 
lens imaging configurations, highlighting the role played 
by phase-sensitive coherence. 

II. SECOND-ORDER COHERENCE 
PROPAGATION 

The imaging configurations we shall consider later ac- 
quire the far-field diffraction pattern or the transverse im- 
age of a transmission mask placed at the source's output 
plane. In both cases linear optical elements — a transmis- 
sion mask, a lens, and a polarizing beam splitter — plus 
free-space propagation lie between the source and image- 
acquisition planes. Because the images are acquired via 
photocurrent correlation, their properties are determined 
by fourth-order correlations of the detected fields. To 
determine these correlation functions, it is necessary to 
propagate the fourth-order correlation function from the 
source plane to the appropriate image-acquisition planes. 
Fortunately, zero-mean Gaussian source states are com- 
pletely determined by their phase-insensitive and phase- 
sensitive, second-order auto- and cross-correlation func- 
tions. Therefore, we need only consider second-order co- 
herence transfer. 

Let t)e~'"''* denote a scalar, z-propagating, 

positive-frequency field operator with center-frequency 
Wo, and ^photons/m^s units. The commutators for the 
baseband envelope at a fixed transverse plane are given 

by El 

[4(Pl,tl),4(P2,<2)]=0 (1) 

ii), El{p2. h)] = S{p^ - p,)d{t2 - ii) , (2) 

where d{-) is the impulse function. 

Free-space paraxial propagation is governed by the 
Huygens-Fresnel principle 12], which states that the 
baseband field operator at z = i is related to the field 
operator at z = by the following superposition integral: 




X hLip- p',LUo + n), (3) 

where 

f4p,0) = j dt E,{p,t)e'^' (4) 

is the Fourier transform of the baseband-envelope field 
operator and 



is the Huygens-Fresnel Green's function for paraxial 
diffraction at frequency w with c being the vacuum light 
speed. 

For Eqs. ©-© to be consistent with Eq. (g]), the in- 
tegrals in Eqs. ([3]) and ([4]) must all be over infinite lim- 
its. This poses no difficulty for the spatial integration, 
which can be taken over the entire z — Q plane, but there 
is a problem with the frequency integration. Because 
Ez{p,t) is a positive-frequency field operator, the lower 
limit of integration in Eq. ([3]) should be — wq- But, if 
we impose this propagation condition we do not preserve 
the commutator in Eq. i.e., assuming that this com- 
mutator applies at z = and employing the diffraction 
integral does not recover the delta-function commutator 
at z = L. In almost all quantum optics situations — 
both theoretical and experimental — we can circumvent 
this issue as follows. If the source at z = only excites 
frequencies that are within some bandwidth ±i7o about 
luq and the measurements performed at z = L are only 
sensitive to frequencies within that bandwidth, then so 
long as r^o < Wo, we can allow the frequency lower limit 
in Eq. ^ to be — oo. In the quasimonochromatic cases to 
be considered below the preceding condition will be sat- 
isfied. However, we will impose the finite lower limit on 
the frequency integral when we address thin-lens imaging 
at the ultimate theoretical limit of broadband SPDC. In 
this case, the relevant commutators are 

[£;(Pi,r!i),£,(p2,^^2)]-0 (6) 
[4(Pi, r!i),4t(p2, ^2)] = 2^5{p2 - p,)S{^2 - ^i) ■ (7) 

The non-Hermitian baseband field operator (p, t) 
has two second-order correlation functions, namely the 
(normally-ordered) phase- insensitive correlation function 

Ki^\p,,P2,t2-h) = {El{p„h)E,{p2,t2)), (8) 
and the phase-sensitive correlation function 

X(f)(Pl,P2,t2-il)= (^.(Pl,tl)^.(p2,t2)), (9) 

in which we have assumed, for simplicity, that the base- 
band field operator is in a complex-stationary state, viz., 
the correlation functions depend on the time difference 
t2 — ti, but not on the absolute times. 

We will find it convenient and insightful to work with 
the frequency spectra associated with Eqs. ^ and ([9]), 
defined as the Fourier transforms, 

Si^\p,,P2,n)^ JdTKi^Hp„P2,T)e^''\ (10) 

for X — n,p. The phase-insensitive and phase-sensitive 
correlation spectra at z = L can be expressed in terms of 
the correlation spectra at z = by evaluating Eq. ITUl) for 
the propagated field operators, via Eq. ([3]), which yields 
the following phase-insensitive spectrum, 

si-\p„P2,n) = II dp[ dp'2 (p'l ,p'2,n) 

X hlipi- p[,UJo + fl)hL{p2- P2,^0 + ^) , (11) 
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and the following phase-sensitive spectrum, 

(Pi , P2 , 17) = J J dp[ dp', s^J"^ ip[ ,p',,n) 

X - Pi,t^O - f7)/lL(P2 - P2:^0 + , (12) 

where * in Eq. (jlip denotes complex conjugation. Note 
that the phase-insensitive correlation spectrum is a 
monochromatic equation, i.e., the frequency dependence 
is + on both sides of the equality, whereas the 
phase-sensitive spectrum is a bichromatic equation in- 
volving ujq ±n. This difference occurs because complex- 
stationary phase-insensitive correlations in time have 
uncorrelated frequency components, whereas complex- 
stationary phase-sensitive correlation functions in time 
have nonzero (phase-sensitive) correlations between fre- 
quency components with equal and opposite detunings 
from the field's center frequency p^ . 

Consider a quasimonochromatic state in which the pre- 
ceding spectra are only nonzero for |il|/(jJo 1- It follows 
that the Huygens-Fresncl principle simplifies to 



El{p, t)= J dp' En{p\t- L/c)hL{p- p',ujn), 



(13) 



and the propagation Green's function in Eqs. (jlip and 
(fT^ become /il(p,u;o ± fi) « /ii(p,wo)- Let us now 
review the far-field propagation regime for this quasi- 
monochromatic situation. We will assume that the field 
at the z = plane is in a zero-mean state with cross- 
spectrally pure, Schell-model correlation spectra [l3| 
given by 



st\p,,p,,n)^T\p,)T{p,)G^"\p,-p,)s^"\n), 



(14) 



S^P\p„p„n)^T{p,)T{p,)G(p\p,-p,)S^P\n). (15) 

With no loss of generality, we require |r(p)| < 1, so 
that it may be regarded as a (possibly complex-valued) 
spatial attenuation of an optical field operator in a ho- 
mogenous and stationary state with separable phase- 
insensitive spectrum S'(")(fi)G^"^(p2 — Pi) and phase- 
sensitive spectrum S''^^-' (r2)G'^P-' (p2 ~ Pi)- This spatial 
attenuation will become the transmission mask to be im- 
aged when we turn our attention to the Fourier-plane and 
thin-lens imaging configurations. For now, however, it is 
convenient to lump this mask together with the source. 

Our primary interest is in sources with narrow G^^"^ (p) , 
for X — n,p, such that T(p) does not vary appreciably 
within a (phase-insensitive or phase-sensitive) coherence 
area. For this case, we may approximate the source cor- 
relation spectra as follows: 

5(")(pi,P2,17) « |r(pjf G(")(p,)5(")(l]), (16) 

si,'\p„p„n) ^T\p^)G^pHp,)s^pHn) , (i?) 



in terms of the sum coordinate p^ = {p, + Pi)/2 and 
the difference coordinate p^ = p, — Pi- This approxi- 
mation simplifies the subsequent analytic treatment con- 
siderably, without significant impact on the fundamental 
physics. The z = L spectra then become 



st\p,.P2,^) 



c^g5(")(ri) 
(27rcL)2 



dp's J dp'd e 



x|r(p^)pG(")(p:,), (18) 



and 



^"^'"'(n n m - ^^i^^!Mp-o(2L^+|pJ^+|p,|V4)/cL 



/ p-i'^o(2p^-pl,+Pd-Pd/2)/ci 



^ e^"«(l''^l'+l''^l'/4)AiT2(p^)G(P'(p:j), (19) 

respectively, where the quasimonochromatic assump- 
tion has permitted our approximating the frequency- 
dependent leading coefficients by their values at the cen- 
ter frequency (l5| . 

Let po denote the coherence radius of the source, which 
we shall assume to be the same for both phase-insensitive 
and phase-sensitive correlations, i.e., po is the radius 
within which G^^^(p), for x — n,p, differ appreciably 
from zero. Also, let oq S> Pq denote the transverse ra- 
dius of |T(p)|2, which will be the photon- flux density 
radius of the source's state just after the transmission 
mask in our imaging configurations. In far field phase- 
insensitive correlation propagation, which applies when 
^oao/5o/2cL ^ 1, the phase term e^'^"^''^'^/'^^ can be ne- 
glected in Eq. (fTS]). and we find that 



4"'(Pl,P2,l]) = 



c^g5(")(r!) 



(27rcL)2 

where T„(A;) and Q^'^'^{k) are the 2-D Fourier transforms, 



r„(fc)^ Jdp'e-^'' '''\T{p')\\ 
g^"\k) = J dp'e-'^ p'G^''\p') . 



(21) 



(22) 



The Fourier-transform duality between the source-plane 
and the far-field phase-insensitive correlation spectra — 
seen in Eq. (PO)) — is the well known van Cittert-Zernike 
theorem for phase-insensitive correlation propagation 
[l^ . A similar duality is present between the source- 
plane and the far-field phase-sensitive correlation spectra, 
but the far-field regime — in which the quadratic phase 
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terms of the integrand in Eq. (fT9|) become neghgible — 
corresponds to ujQa'^l2cL <C 1, which is more stringent 
than the far-field condition for the phase-insensitive case. 
In this regime, we find that 



■ L^m free space propagation (far fiel^) 



(27rcL)2 



(23) 



gives the far-field phase-sensitive correlation spectrum, 
with 



and 



Tp(fe) = jdp'e-'^p'T\p') 
g^P^k) = J dp'e-'^ p' G^P\p') . 



(24) 



(25) 



By analogy to the phase-insensitive case, we refer to the 
Fourier transform relation in Eq. ((23|) as the van Cittert- 
Zernike theorem for phase-sensitive coherence propaga- 
tion. 

To conclude our review of far-field coherence prop- 
agation it is worth emphasizing the similarities and 
differences between Eqs. (P0|) and (^5]) . The source- 
plane transmission mask, T{p), has been assumed to 
be a slowly-varying and broad function in comparison 
to the rapidly decaying G'^^\p). Thus, for x = n,p, 
Fourier-transform duality implies that Tx{k), decays 
more rapidly than fj^^^ (fe). Therefore, the far-field phase- 
sensitive correlation function consists of a narrow func- 
tion of Pg multiplying a broad function of p^, whereas 
the corresponding phase-insensitive correlation function 
consists of a narrow function of p^ times a broad func- 
tion of Pg. Owing to this difference, point pairs on the 
transverse plane in the far field that are symmetrically 
disposed about the origin, viz., points satisfying « 0, 
have appreciable phase-sensitive correlation. The phase- 
insensitive correlation, however, is highest between point 
pairs that are in close proximity on the transverse plane, 
i.e., point pairs obeying « 0. In addition, if we eval- 
uate the correlations at a single transverse point, i.e., 
when p^ = 0, we find that the phase-insensitive cor- 
relation traces out the broad envelope Q'^^^ujqPs/cL), 
whereas the phase-sensitive correlation traces out the 
narrow function 7p{2ujoPg/cL), a property we shall make 
use of in the following section. Finally, it is relevant to 
emphasize that Eq. (^5)1 is a general property of phase- 
sensitive coherence propagation that applies regardless 
of whether the source state is classical or nonclassical 



III. TWO-PHOTON DIFFRACTION PATTERN 
IMAGING 

Consider the experimental setup, shown in Fig. [1] 
whose purpose is to obtain the far-field diffraction pat- 
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FIG. 1: (Color online) Imaging the far-field diffraction pattern 
of a transmission mask. PBS, polarizing beam splitter. 



tern of a source-plane transmission mask of (possibly 
complex- valued) field transmissivity T{p). An exper- 
iment using this setup — which we will explain in de- 
tail shortly — was reported in [3| as a proof-of-principle 
demonstration of biphoton-state quantum lithography. 
That experiment exhibited a factor-of-two compression 
in the fringe pattern produced using a two-slit transmis- 
sion mask as compared to what was obtained with con- 
ventional (coherent-state) illumination at the same wave- 
length of that two-slit mask. The observed fringe-pattern 
compression was therefore interpreted as the expected 
result for quantum lithography using a pair of entan- 
gled photons Our aim in this section is to show 
that it is phase-sensitive coherence, not entanglement 
per se, that is responsible for this fringe-pattern com- 
pression. In particular, we will find that classical phase- 
sensitive light and biphoton-state light yield identical im- 
ages, except for the image produced with the classical- 
state source being embedded in a prominent featureless 
background that is absent for the case of biphoton illumi- 
nation. Thus, the fringe-pattern compression previously 
ascribed to the entangled nature of the signal and idler 
photons is in fact due to their having a phase-sensitive 
cross correlation, which is a property that classical states 
can also possess. Signal-idler entanglement, which in 
the Gaussian-state framework is a stronger-than-classical 
phase-sensitive cross correlation between the signal and 
idler fields, is responsible for dramatically improving the 
contrast of this image when the illumination is broad- 
band. 

Now let us flesh out the details of the preceding asser- 
tions within the context of the Fig. [1] setup. The source 
will be taken to be frequency-degenerate type-II phase- 
matched SPDC with a continuous-wave pump. It gener- 
ates paraxial, z-propagating signal (S) and idler (/) fields 
in orthogonal polarizations with common center fre- 
quency luq. Their positive- frequency field operators will 
be denoted as Es{p,t)e~'''^°'^ and Ej{p,t)e-'''^°* . With 
a nondepleting plane- wave pump and ignoring boundary 
effects due to the nonlinear crystal's finite cross-section, 
these two output fields are in a zero-mean, jointly Gaus- 
sian state that is homogeneous and stationary. More- 
over, the signal and idler will then have identical fluores- 
cence spectra and maximum phase-sensitive cross corre- 
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lation, but no phase-sensitive autocorrelation or phase- 
insensitive cross correlation [J, It follows that after 
passing through the transmission mask the jointly Gaus- 
sian state is fully determined by the Schell-model auto- 
correlation spectra [14] 



(26) 



for m = S, I, and the phase-sensitive cross-correlation 
spectrum 



4^|(Pl,P2,l])=5r(Pl,P2,f^), 



(27) 



where S*^"^ and S^^^ are given by Eqs. (jl6p and (|17p 
respectively, and we shall assume that they satisfy the 
quasimonochromatic condition. 

In the Fig. [T] setup, the signal and idler fields both 
propagate over an L m free-space path — assumed to 
be sufficiently long to satisfy the far-field condition, 
a;oao/2cL <C 1, for the phase-sensitive source — and then 
are separated by a polarizing beam splitter such that they 
impinge on separate pinhole detectors, each centered on 
the transverse-plane coordinate p with respect to its op- 
tical axis. Because linear transformations of zero- mean 
Gaussian states are still zero-mean Gaussian [l| , and be- 
cause free-space diffraction is a linear transformation, we 
need only determine the second-order moments at the 
detection planes to determine the joint state of Ei{p,t) 
and E2 {p, t) , which denote the far- field propagated field 
operators of the signal and the idler respectively. Thus, 
quasimonochromatic paraxial diffraction into the far field 
results in phase-insensitive autocorrelation spectra given 
by Eq. ([20)1 and a phase-sensitive cross-correlation spec- 
trum given by Eq. ^ 

The two pinhole photodetectors are assumed to have 
identical parameters: quantum efficiency 77; photosensi- 
tive area A; and current-pulse output qhsit) from detec- 
tion of a single photon, where q is the electron charge 
and JdthB{t) = 1. Then, the time-average photocurrent 
cross-correlation at the detection planes has an ensemble 
average (9l. il9ii] 



C{p) 



T 



T/2 



T/2 



(28) 



in which the equal-time photocurrent cross correlation is 
given by 



(*i(i)«2(<)) = g J J duidu2 

X {e\{p, ui)eI{p, U2)Ei{p, ui)E2{p, U2)) 

X heit - ui)hBit - U2) . (29) 

Here, we have approximated the integrals over the pin- 
hole detectors' photosensitive regions by the value of the 
integrand at p times A^. For our Gaussian-state Ei and 
E2, the fourth-order field moment in Eq. (^5)1 reduces 



to a sum of products of second-order correlation func- 
tions by virtue of the Gaussian moment-factoring theo- 
rem [3, . This procedure simplifies the photocurrent 
cross-correlation expression to 



Cip) = Coip) + Cp 



/ 2ljop\ 



(30) 



Here, 



Co(p) 



47r2c2i2 



r„(o)(?(")f^^ 



V cL 



hsit) 



(31) 



is a non-image-bearing background, which is broad and 
featureless owing to Q'^"'\k) being the Fourier trans- 
form of the narrow spatial-domain correlation function 
G(")(P)- The second ter m in Eq. (|30p is the image- 
bearing term. Its constant factor is 



^ 14^2^27^2 



\T-^{S^P\n)}\^i.hBi^'hB 



(32) 

where J- ^{-j denotes the inverse Fourier transform of 
the bracketed term [see Eq. ([TOl) for our Fourier transform 

sign convention], ★ denotes convolution, and hg repre- 
sents the time- reversed impulse response. From Eq. (j30l) . 
we see that the image-bearing term is proportional to 
\Tp{2ujop/cL)\'^, which is the far-field diffraction pattern 
of the square of the mask's field transmissivity T{p). 

Let us compare the imaging characteristics of this im- 
ager to those of a conventional classical imager that 
utilizes a coherent-state beam to illuminate the mask, 
and a single (scanning) pinhole detector located in the 
far field that records the diffraction pattern. If we as- 
sume the field impinging on the transmission mask is a 
monochromatic plane wave at center frequency t^o a-nd 
with photon-flux density /q, the field just after the mask 
is in the coherent state satisfying 



(33) 



Ep{p,t)WhT{p)) - ^hT{p)WhT{p)) 



Because free-space propagation is a multimode beam 
splitter relation ;12], the detection-plane field operator 
Ei{p,t) is also in a coherent state, whose eigenfunction 
is determined by substituting ^/T(jT{p) into the classical 
Huygens-Fresnel diffraction integral, i.e., Eq. ^ with the 
field operator replaced by the coherent-state eigenfunc- 
tion. 

We shall assume the path length L satisfies the far-field 
condition, ujoaQ/2cL ^ 1, for coherent-state diffraction 
[10], so that the quadratic phase term in the Huygens- 
Fresnel Green's function becomes negligible and the mean 
photocurrent becomes 



m) 



47r2c2L2 



/o / dthB{t) 



(34) 
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where there is no background and the image term is given 
by, 



(35) 



If T{p) only takes values zero or one — as was the case 
for the two-slit transmission mask employed in Q — then 
r^(p) — T{p) and the biphoton source yields a far- 
field diffraction pattern proportional to |7^(2a;op/cL)p, 
whereas the coherent diffraction pattern is proportional 
to \Tp{u!op/cL)\'^. Thus, the far-field pattern observed 
with the biphoton source is spatially compressed by a 
factor of two relative to that obtained when the coherent- 
state source is employed, which is why the biphoton case 
has been said to beat the classical resolution limit. How- 
ever, it is worth re-emphasizing that Eq. (j30[) is true for 
both classical and quantum Gaussian-state sources whose 
signal and idler fields have a nonzero phase-sensitive cross 
correlation. Consequently, it is the phase-sensitive co- 
herence of the source — and not the entanglement of the 
biphoton — that is the fundamental cause for the factor- 
of-two compression in the far-field diffraction pattern 
as compared to what is obtained using a coherent-state 
plane wave. Finally, it is important to note that unless 
T^{p) oc T{p), as we have above when T{p) is a zero- 
one function, then the diffraction pattern acquired from 
phase-sensitive sources will be distorted relative to what 
will be obtained with a coherent-state field. 

Utilizing phase-insensitive Gaussian-state light in the 
Fig. [T] configuration does not result in a photocurrent 
cross correlation containing a diffraction-pattern image. 
This imaging failure occurs because, from Eq. (PO)) . the 
equal-position correlation in the Fourier plane traces out 
Qq"'\ujq p / cL) , which does not contain any information 
about the transmission mask T{p). However, by modify- 
ing the Fig. [1] setup to have one detector scan —p, while 
the other scans p, then the photocurrent cross correlation 
will contain an image of \Tn{2uJap/cL)\'^. For a zero-one 
function transmission mask, this phase- insensitive imager 
also achieves the factor-of-two pattern compression in 
comparison with the image formed using a coherent-state 
plane wave. In this regard we note that phase-insensitive 
Gaussian-state sources are classical Q . 

Thus far we have considered only the image-bearing 
term in Eq. ([50)1 . Now we will address the image contrast. 
For simplicity, we will assume that T(p) is real valued. 
In addition, we restrict ourselves to an observation region 
TZ that encompasses the image-bearing term in Eq. (I30p . 
and we define the contrast as 



C 



max7^[C(p)] - minu[C{p)] 
Co(0) 



(36) 



so that the numerator yields the dynamic range of the 
image-bearing term in the photocurrent correlation C (p) , 
while the denominator is the featureless background. 

Here we compare the contrast from classical and quan- 
tum sources that have identical autocorrelation spec- 
tra and the maximum phase-sensitive cross correlation 



allowed in classical and quantum physics, respectively. 
When the source is in a classical Gaussian state, whose 
autocorrelation spectrum — ^just before the transmission 
mask— is g'-''\k)S'-''^n), the maximum magnitude for 
the phase-sensitive cross-correlation spectrum is equal to 
the autocorrelation spectrum [1], i.e., the Gaussian state 
with maximum classical phase-sensitive cross correlation 
satisfies 



(37) 



Taking the phase of this phase-sensitive spectrum to be 
zero, and recalling that J dthB{t) = I, the contrast with 
classical phase-sensitive Gaussian-state sources can be 
written in the form 



where the spatial (s) factor is given by 

maxfc[|rp(fc)|2] - minfc[|rp(fc)|2] 



with equality if T{p) is real, so that 



< I , 



(38) 



(39) 



,(c) 



|j^-H5(")(f7)}|^*/iB* he 



(/dOS'(")(f7)/27r 



(40) 



for such masks. For analytical convenience, let us take 
the spectral part of the phase-insensitive autocorrelation 
function to be Gaussian with e~^-attenuation baseband 
bandwidth 2/To, i.e.. 



(41) 



and let us take the baseband impulse response hB{t) to 
be a Gaussian with e"^-attenuation time duration T^, 



VIZ., 



hsit) 



(42) 



With these assumptions, we find that the classical con- 
trast for a real- valued mask is 



I 



v/1 + (7d/2To)2 



(43) 



which is approximately unity for narrowband sources 
that satisfy Td ^ Tq. On the other hand, in the broad- 
band limit Td ^ Tq, we have 



C(^) « 2To/Td « 1 



(44) 



so the contrast is severely degraded in this case. 

Now consider a nonclassical Gaussian state with the 
maximum phase-sensitive cross correlation. In the low- 
brightness regime, i.e., when 5'(")(r2)e;^"^(fe) < 1, the 
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maximum phase-sensitive cross-correlation spectrum is 
approximately ^] 

\s^p\n)g^p\k)\ » (45) 

which is much higher, in this limit, than the classical 
maximum given by Eq. (j37p . Taking the phase of this 
correlation to be zero, the contrast is found to factor into 
the product of spatial and temporal terms, with ci*'' = 

(c) 

Cs , and the temporal term given by 



-><- 



.(■7) 



^(")(0) (/dr2S'(")(17)/27r) 



(46) 



Once again using Eq. (j4ip for the fluorescence spectrum 
and Eq. (|42|l for the baseband current filter, we obtain 



Cf") 2/g(")(0)5(")(0) V 1 + r|/2r2 (47) 

for real- valued T{p). Here, the narrowband contrast 
C(i) = 2/e;(")(0)S'(")(0) is very high because of the low- 
brightness condition, and even for broadband fields the 
contrast. 



c(«) = 2V2ro/rdg^("no)5^"^(o) 



(48) 



may be high. In particular, in the biphoton regime, 
wherein t/'-"'' (0)Td <C 1 (low flux as well as low bright- 
ness) also prevails, very high contrast is predicted in 
this broadband limit [2l!|, which is in agreement with 
the background- free diffraction pattern reported in Q. 
Therefore, low-brightness quantum Gaussian-state fields 
have a contrast advantage over classical phase-sensitive 
Gaussian-state fields when the phase-sensitive cross cor- 
relation is measured via a photocurrent correlation mea- 
surement, and the biphoton state yields images with neg- 
ligible background even when it is a broadband state. 

In summary, in this section we have studied a corner- 
stone proof-of-principle experiment for quantum optical 
lithography by applying Gaussian-state analysis and the 
coherence theory results from Sec. |lT] to the propaga- 
tion of classical and quantum phase-sensitive cross cor- 
relations. Our analysis has shown that the only perfor- 
mance difference between using a biphoton-state source 
and a classical phase-sensitive Gaussian-state source in 
the Fig. 1 setup is the diffraction-pattern image's con- 
trast. The resolution improvement seen with a bipho- 
ton source is entirely due to the diffraction properties 
of the phase-sensitive cross correlation between the sig- 
nal and idler fields, hence it is also achievable with a 
pair of classical Gaussian-state fields with phase-sensitive 
cross correlation. However, low-brightness quantum 
sources achieve higher contrast than classical sources, 
which permits imaging with broader bandwidth quan- 
tum sources. Finally, the broadband biphoton state 
yields very high contrast images, which is the reason why 
biphoton-state quantum lithography experiments have 
yielded background-free diffraction-pattern images Q. 
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FIG. 2: (Color online) Two-photon thin-lens imaging of a 
transmission mask. PBS, polarizing beam splitter. 



IV. TWO-PHOTON THIN-LENS IMAGING 

Let us now consider using an optical source with low 
spatial coherence and a thin lens to image a transmission 
mask placed at the source plane, as depicted in Fig. ^ 
Primary attention in our analysis of this experimental 
setup will be given to the resolution limitations imposed 
by the finite aperture of the lens, as previous work has 
claimed that a factor-or-two resolution improvement ac- 
crues when a broadband biphoton source is employed 

As in the previous section, we shall assume an SPDC 
source that generates zero-mean Gaussian-state signal 
and idler beams whose phase-insensitive correlation spec- 
tra at the exit plane of the transmission mask are given 
by Eq. P5|) and whose phase-sensitive cross-correlation 
spectrum at that plane is given by Eq. (|77|) . The opti- 
cal fields first propagate through a di-m-long free-space 
path according to Eq. ([3]). A thin lens of radius R and 
focal length / embedded in an otherwise opaque screen is 
placed on this plane. This thin lens will be assumed to be 
entirely free of chromatic aberration over the frequency 
range of interest, i.e., each frequency component of the 
impinging field is multiplied by circ(|p|/i?)e~''^l''l ^■^'^■f ^ 
where uj is a passband frequency centered around ujq and 
the circle function is 



circ(a;) 



1, k|<l, 
0, otherwise. 



(49) 



Finally the field at the exit plane of the lens propagates 
d2 m in free space to reach the image plane, which is 
defined by the lens maker's formula -I- = f~^. 
The image-plane signal and idler fields are separated by 
a polarizing beam splitter, after which each illuminates a 
pinhole photodetector located at transverse coordinate p 
relative to its optical axis. The resulting photocurrents 
are then correlated to obtain the same fourth-order field 
measurement given in Eqs. ([^5]) and in terms of the 
image-plane field operators Ei{p,t) and E2{p,t). 

The overall mapping from the source-plane field op- 
erators to the image-plane field operators is linear, and 
therefore Ei (p, t) and i?2(p, t) are in a zero-mean jointly- 
Gaussian state. With the simplifying assumption that 
the detectors have electrical bandwidths much broader 
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than the source spectra [2^ , we can combine and reduce 
Eqs. (pS)) and for the biphoton state to obtain 

C{p) « <z\2^2 ^^t i)^t ^)^^ i)^^ (50) 



<2^(P,P,0)|^ 



(51) 



where 



kI;:\{p,P,t) ^ {El{p,t)E,{p,t + r)), (52) 



P, r) = (^„(p, t)i?,(p, t + r)), 



(53) 



for m,£ = 1,2, and Eq. (1511) follows from the Gaussian 
moment-factoring theorem ^16i. ,19.]. Furthermore, as we 
have determined in the previous section, for maximally- 
entangled Gaussian states with low-brightness and low- 
flux, i.e., the biphoton state, the second term in (|5ip is 
much stronger than the first, permitting the approxima- 
tion 



Cip) 



l!'i(p,P,o)p. 



(54) 



Therefore, for biphotons a photocurrent correlation mea- 
surement with broadband detectors is a means for mea- 
suring the squared magnitude of the phase-sensitive 
cross correlation between the image-plane field operators 
p^ . In this section we demonstrate that the interesting 
biphoton-state results predicted for this imaging configu- 
ration are a consequence of the phase-sensitive cross cor- 
relation, and the photocurrent correlation does not play 
a role beyond facilitating its measurement. Hence, in 
the remainder of this section we shall bypass this pho- 
tocurrent correlation measurement and focus directly on 
the phase-sensitive and phase-insensitive cross correla- 
tions between the image-plane field operators. 

The frequency-domain image-plane field operators are 
given by the following linear transformation of the 
frequency-domain source-plane field operators. 



£m{p, ^) y '^P' ^(P^ P' ^ ^0 + '^)£l{p' , ^) + Cm{P, ^) , 

(55) 

for (£, to) = (S*, 1), (/, 2), where £m{p,^) is an auxil- 
iary vacuum-state operator such that Em{p,t) satisfies 
the free- field commutators given in Eqs. ([T]) and 
The point-spread function /i(p, p',u;), found from the 
Huygens-Fresnel principle and the lens transfer function, 
is given by 



h{p,p\iu)^n(r{p,p'),u;/u;o)e^^^P-P''-\ (56) 



where 



and 



r(p,p')^^^MipM + p'| 



(57) 



(58) 



with 2Ji{x)/x for a; > being the well-known Airy func- 
tion. The phase term in Eq. ([55]) is 

(j){p, p', L0)^L0 {dl +d2 + \pf/2d2 + |pf /2di) /c, (59) 

which incorporates the group delay arising from the 
{dl + d2)-m propagation, and the parabolic phases at 
the source and image planes that are associated with the 
diffraction process. 

Because our concern is with the resolution limit im- 
posed by the lens having a finite radius i?, we will further 
simplify our analysis by assuming spatially-incoherent 
source statistics and appropriate focusing at the source 
plane to compensate for the parabolic phase in Eq. (|59p . 
These assumptions simplify the phase-insensitive au- 
tocorrelation functions and the phase-sensitive cross- 
correlation function of the /q photons/m^s signal and 
idler fields — given in Eqs. ITBl) and (|17p — to 

5(")(Pi,P2,17) = |T(pi)|2[27rc/(c^o + r!)]'/o 

X 5(p2-pi)s(")(f}) (60) 

and 

4'')(Pi,P2,l]) = e— "'"^'^/^''^ T\p,) 

X [{2ncr/{ojl-n^)]l,6{p,-pi)s('^\n), (61) 

respectively, where s'-^^(r2)/27r, for x — n,p, are nor- 
malized (unity area) spectra |23 . Evaluating the phase- 
insensitive autocorrelations and the phase-sensitive cross 
correlation of the two image-plane fields at equal spatial 
coordinates (relative to their optical axes) and at equal 
times, yields 

K(^l{p,p,0)^ f dp'|T(p')|V(KP,P')), (62) 

for m = 1, 2, and 

k[^1{p,P,0) = e*"o(2di-h2<i2-H|p|Vd2)/c 

X / dp'T\p')gp{rip,p')). (63) 



The point-spread function in the superposition integral 
involving |r(p)p is 

X \n{r,l + n/LUo)\yil + n/LUo)\ (64) 

and it determines the phase-insensitive autocorrelation 
functions. Likewise, the point-spread function in the su- 
perposition integral involving r^(p) is ^] 

X H{r, 1 + n/LUo)n{r, 1 - n/ujQ)/{l - n^/LU^),{65) 
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and it controls the phase-sensitive cross-correlation 
function. Therefore, apart from an unimportant 
parabolic phase factor, the most important difference 
between phase-insensitive and phase-sensitive coherence 
propagation — insofar as two-photon thin-lens imaging is 
concerned — is the frequency coupling between ztfi/wo 
that is present in Eq. ([65|l but is absent from Eq. ([64|) . 
In the quasimonochromatic limit, however, this coupling 
becomes insignificant, because 



(66) 



so that 
s(p)(f7). 



gnir) — gp{r) prevails whenever s'^"^(il) — 
Thus the quasimonochromatic point-spread 
function for the phase-insensitive correlation is identical 
to its quasimonochromatic phase-sensitive counterpart. 

However, 5„(r) and gp{r) begin to differ as the band- 
width of the source increases. Suppose that the normal- 
ized phase-insensitive and phase-sensitive source spectra 
are both taken to be flat over a 2M^-bandwidth window. 




\n\ < w, 

otherwise. 



(67) 



Substituting this expression into Eqs. ([64]) and ((65)) per- 
mits us to express the point-spread functions as the fol- 
lowing dimensionless integrals. 



gnir) 



' du ■ 



2c^dldlW 



and 



W/uio 



2c^dldlW 



(68) 



gpir) 



W/ujo 



du 



W/uo 

Ji{r{l + u)) Ji{r{l~u)) 



(69) 



In the quasimonochromatic limit, in which W/loq 1 
holds, both point-spread functions simplify to 



,{r) = gp(r) 



IpL^^R* ^2Ji(r)^2 
Ac^djdj 



(70) 



Hence with a quasimonochromatic source there is no dif- 
ference between the image of a real-valued transmission 
mask acquired with phase-insensitive (thermal) illumi- 
nation or phase-sensitive (classical or quantum) illumi- 
nation. 

The point-spread functions for broader bandwidth 
sources are plotted in Figs. [3] and [4] at two different W 
values. The W = ujq / 4: phase-sensitive point-spread func- 
tion, shown in Fig. [31 represents imaging performance for 
unusually broadband SPDC 25|. The W ~ ujo phase- 



sensitive point-spread function, plotted in Fig. [H rep- 
resents imaging performance at the ultimate theoretical 
limit of SPDC bandwidth. The point-spread functions in 
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FIG. 3; (Color online) Comparison of the imaging point- 
spread functions for phase-insensitive (PIS) and phase- 
sensitive (PS) correlations when W/loo = 0.25, and when the 
imaging source is quasimonochromatic (QM). The normaliz- 
ing coefficient is «: = /oi^o^^/4c^rfirf2- 



these figures show that the peak amplitude of the phase- 
insensitive function increases to 



'b 



(71) 



whereas that of the phase-sensitive point-spread function 
attenuates to 



5f>(0) 



Ac^djd 



1 - 



3^0 



(72) 



relative to the peak amplitude in the quasimonochro- 
matic limit as the source bandwidth increases. The 
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FIG. 4: (Color online) Comparison of the imaging point- 
spread functions for phase-insensitive (PIS) and phase- 
sensitive (PS) correlations at the ultimate theoretical limit 
of SPDC bandwidth (W/t^o = 1), and the quasimonochro- 
matic (QM) limit. The normalizing coefficient is k = 



(1 -|- uf' factor multiplying the frequency-resolved Airy 
patterns in Eq. where |u| < Wj^Q < 1, is responsi- 
ble for the increase in (0) with increasing source band- 
width. This scaling favors the 6/ue-detuned frequency 
contributions to the phase-insensitive point-spread func- 
tion. So, because of the quadratic scaling, the average 
of the amplitude increase for m > versus the ampli- 
tude decrease for u < is greater than one. Thus, the 
peak value of gn(0) increases with increasing source band- 
width. On the other hand, the scaling for the frequency- 
resolved Airy patterns in the integrand of Eq. ([M)) is 
(1 — w)(l + u) = 1 — M^, implying that all detuned fre- 
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red- 
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FIG. 5: (Color online) Image-plane spot diameters for differ- 
ent frequency components of a broadband point source. 



quencies are attenuated, which results in .gp(0) decreasing 
as the source bandwidth increases. 

Figures El^a) and SJa) indicate that some narrowing 
of the main lobe occurs for both the broadband phase- 
sensitive and phase-insensitive point-spread functions, 
with the latter suffering from a slower-decaying tail. 
From a practical perspective this main-lobe narrowing 
behavior is of little interest. Taking the resolution to be 
set by the first zero in the point-spread function, Fig.IS^b) 
shows that the resolution benefit offered by broad- 
band phase-sensitive as compared to quasimonochro- 
matic (phase-sensitive or phase-insensitive) imaging is 
merely a factor of 1.14. From this figure we also note 
that the tail of the phase-insensitive point-spread func- 
tion traces the envelope of the oscillations in the phase- 
sensitive point-spread function, so no appreciable loss 
of resolution results from the former's slowly-decaying 
tail. Even at the W = loq ultimate bandwidth limit, 
the resolution improvement offered by the broadband 
phase-sensitive point-spread function — as found from 
Fig. m^b) — is only a factor of 1.38. In this limit, the tail 
of the broadband phase-insensitive point-spread function 
falls off somewhat — but not dramatically — slower than 
the envelope of the phase-sensitive point-spread func- 
tion's oscillations. 

The difference between the behavior of the phase- 
insensitive and phase-sensitive point-spread functions as 
a function of the source bandwidth deserves closer ex- 
amination to understand the underlying physics. Recall 
that a source generating a complex-stationary baseband 
field around a center frequency wo is a superposition 
of monochromatic field components which have phase- 
insensitive autocorrelations at each frequency and phase- 
sensitive cross correlations between frequencies that sum 
to 2loq. Thus, the phase- insensitive correlation, mea- 
sured at a given spatiotemporal coordinate {p,t), is a 
superposition of all the different autocorrelations at de- 
tunings Q over the fluorescence bandwidth of the source. 
On the other hand, the phase-sensitive correlation mea- 
sured at (/o, t) is a superposition of all the cross corre- 
lations between frequency components detuned by ifi, 
over the phase-sensitive bandwidth of the source. Now, 
consider a point source at the source-plane that emits sig- 
nal and idler fields that have nonzero phase-insensitive 
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autocorrelations and phase-sensitive cross correlation. 
From Eq. (|55p . signal and idler frequency components 
at uj = luq + will yield image-plane spots with common 
radius cd2/{ujo + ^)R- Thus, as shown in Fig.O lower fre- 
quency components produce broader spots on the image 
plane than do higher frequency components. The phase- 
insensitive autocorrelation — of either image-plane field — 
measured by scanning a point detector on the transverse 
plane, therefore decays slowly as \p\ increases, because 
of the large spots from the lower frequencies. However, 
this slowly-decaying tail does not cause a significant in- 
crease in the point-spread function's width, because the 
quadratic weighting coefficient in Eq. ((68)) accentuates 
blue-detuned frequencies and attenuates those that are 
red detuned. On the other hand, if we are measuring 
the phase-sensitive cross correlation between the signal 
and idler fields (26j . we are in effect measuring the su- 
perposition of the cross correlations between the ujq + 
signal-field component and the ujq — Q idler-field compo- 
nent, where H E [— ljoi'^o]- For > 0, the former yields 
a narrow spot of radius 0^2/ (wq + and the latter 

yields a broad spot of radius cc?2/(wo — ^)R- Because the 
phase-sensitive cross correlation is given by their prod- 
uct, however, the narrower radius from the higher fre- 
quency determines the radius within which there is ap- 
preciable phase-sensitive coherence. Furthermore, this 
coherence radius is symmetric in fl, so, as the phase- 
sensitive bandwidth of the source increases, the width 
of the image-plane phase-sensitive point-spread function 
decreases. However, the weighting coefficient 1 — in 
Eq. ([M)) counteracts this advantage by attenuating the 
frequency contributions with higher detunings, such that 
the net reduction in the main lobe's width is very small. 



Notice that we have made no reference to the classical 
or quantum nature of the source in explaining the physics 
governing the point-spread functions' width. Thus, this 
effect is entirely a consequence of phase-sensitive versus 
phase-insensitive source correlations and scalar paraxial 
diffraction theory, both of which are valid in the classical 
and quantum theories of light. The quantum nature of 
the fields, therefore, does not play a role in determining 
the resolution capabilities of thin-lens correlation imag- 
ing, regardless of whether the source has phase-sensitive 
or phase-insensitive coherence. However, for particular 
measurement schemes, nonclassical field states may offer 
contrast advantages akin to those found in the previous 
section for diffraction-pattern imaging. In particular, in 
this section we have determined that the phase-sensitive 
correlation differs from its phase-insensitive counterpart 
only in the broadband limit. Thus, if we opt to utilize a 
photocurrent correlation measurement, then the contrast 
will be significantly better when the broadband fields' 
state is maximally-entangled (nonclassical) and has low- 
brightness, which encompasses the biphoton state. 



V. DISCUSSION 

SPDC with vacuum-state inputs generates signal and 
idler fields in a zero-mean jointly-Gaussian state, with 
nonzero phase-insensitive autocorrelations and a phase- 
sensitive cross correlation that fully determine their joint 
state. When the output state is driven to the low- 
brightness, low-flux limit, this Gaussian state becomes 
equivalent to a dominant vacuum state plus a weak 
biphoton contribution, in which the biphoton wave func- 
tion equals the phase-sensitive cross correlation between 
the signal and idler fields. On the other hand, clas- 
sical imagers have traditionally utilized optical sources 
in thermal states or coherent states, both of which are 
Gaussian states but have only nonzero phase-insensitive 
correlations. Hence, quantum imaging experiments that 
rely on biphoton sources, as well as conventional classi- 
cal imaging configurations, can be unified and generalized 
by studying the imaging characteristics of Gaussian-state 
sources. Furthermore, such states are fully characterized 
by their first and second moments, and are closed under 
linear transformations on the field operators. So, imag- 
ing configurations utilizing Gaussian-state sources, linear 
optical elements and free-space propagation can be fully 
understood in both the classical and quantum regimes by 
tracking the evolution of the first and second moments of 
the fields from the source plane to the detection planes. 

A particularly relevant distinction that has been over- 
looked in most previous work is the phase-sensitive na- 
ture of the correlation between the two photons in a 
biphoton state, as opposed to the phase-insensitive cor- 
relation that is present between thermal-state fields. 
Phase-sensitive coherence has propagation characteris- 
tics that differ from those of phase-insensitive coher- 
ence. Furthermore, complex-stationary phase-sensitive 
correlations have cross-frequency couplings that are not 
present in complex-stationary phase-insensitive correla- 
tions. Distinctions such as these often underlie the inter- 
esting observations and theoretical predictions in quan- 
tum imaging. However, phase-sensitive coherence is not 
exclusive to nonclassical states (such as the biphoton). 
Classical Gaussian states (random mixtures of coher- 
ent states) may very well have nonzero phase-sensitive 
correlations, and those features in quantum imaging 
that stem from the phase-sensitive coherence between 
the two photons in a biphoton state can be replicated 
with classical phase-sensitive sources, as we have pre- 
viously demonstrated for optical coherence tomography 
and ghost imaging 

In this paper we continued to distinguish the truly 
quantum phenomena in quantum imaging theory and ex- 
periments from the phase-sensitive coherence phenomena 
that can be exploited both in the classical and quan- 
tum regimes. Toward this end, we performed Gaussian- 
state analyses of two significant experimental configura- 
tions in biphoton-state quantum imaging. In Sec. IIIII we 
showed that the factor-of-two spatial compression in the 
far-field diffraction pattern of a transmission mask placed 



12 



at the source plane is precisely due to the phase-sensitive 
cross correlation between the signal and idler fields, in 
both the classical and quantum regimes. Indeed, the 
only significant difference — insofar as this experiment is 
concerned — between phase-sensitive classical and quan- 
tum sources is the image contrast when photocurrent cor- 
relation measurements are employed. Narrowband classi- 
cal Gaussian states can achieve acceptable contrast, but 
the contrast degrades severely when the source is broad- 
band. On the other hand, with low-brightness quan- 
tum Gaussian states that are maximally-entangled, the 
contrast is high for both narrowband and broadband 
sources. Note that the strength of the background in 
the signature may be a relevant factor in determining 
whether a classical or quantum source is more desirable 
for a particular application. For example, in photolitho- 
graphic applications, in which extraneous noise may be 
eliminated by virtue of operation in a controlled envi- 
ronment, a biphoton-state source in combination with a 
two-photon absorber at the detection plane generates an 
optical image with no background, whereas a classical 
phase-sensitive source yields significant background that 
requires postdetection processing prior to etching. Hence 
the contrast advantage offered by the biphoton state — 
which cannot be replicated by classical phase-sensitive 
light — is a desirable feature in this case. 

In Sec. IIVI we compared thin-lens imaging of a 
source plane transmission mask using incoherent phase- 
insensitive light to the same imaging arrangement using 
phase-sensitive light. When the sources are narrowband 
(quasimonochromatic), the point-spread functions of the 
two cases turn out to be identical, yielding no resolution 
difference between the various source possibilities. As the 
source bandwidth increases, the point-spread functions 
for the phase-insensitive and phase-sensitive correlation 
functions become narrower, with the phase- insensitive 
point-spread function developing a more slowly decay- 
ing tail. The differences between the two cases stem 
from complex-stationary source statistics, and frequency- 
dependent free-space diffraction. Once again, the bipho- 
ton state facilitates a high-contrast image and a conve- 
nient measurement apparatus (coincidence counting) for 
detecting phase-sensitive correlation, but it is not respon- 
sible for the physics governing the changes to the point- 
spread functions. 

Although Sec. lIIII concentrated on a biphoton proof-of- 
principle experiment for quantum lithography, the driv- 
ing motivation for quantum optical lithography is the iV- 
fold improvement in etching resolution that is predicted 
for a system using iVOOiV states The A^OOA^ state 

is an equal- weight superposition of two pure states: an 
A^-photon signal field and a vacuum-state idler, plus a 
vacuum-state signal field and an A^-photon idler. The 
NOON state is nonclassical; its P representation in terms 
of coherent states is not a proper probability density. The 
N — 2 case can be achieved with biphoton states, viz., 
the output of a 50-50 beam splitter when the two inputs 
are the signal and idler fields from SPDC operating in 



the low-brightness, low-flux regime. Unfortunately, gen- 
erating NOON states for A^ > 2 has proven challenging. 
Thus far the interference fringes for the A^ = 3 and N = A 
cases have been demonstrated in proof-of-principle exper- 
iments [1^ HO] , showing factors of 3 and 4 fringe compres- 
sion respectively, and efforts to generate higher orders 
continue. NOON states with N > 2 are not Gaussian 
states or any limiting form of Gaussian states, because 
their second-order moments do not determine the state. 
Therefore, the Gaussian-state analysis presented in this 
paper does not generalize to NOON states with N > 2. 
As a result, to better appreciate the fundamental physics 
that leads to improved resolution with these sources, it is 
of great interest to develop a unifying coherence theory 
for higher-order moments of continuous field operators, 
and perform an analysis for these moments to determine 
whether the advantages observed with these states are 
truly due to their nonclassical nature or due to a mea- 
surement of a 2Ath-order moment of the field operator. 

The analyses presented in Sees. IIIII and IIVI reveal 
that the physics governing the resolution improvement in 
Fourier-plane and thin-lens imaging are different. Specif- 
ically, the improvement in resolution that is observed 
in the Fourier-plane measurement is due to the dif- 
ference in paraxial propagation of phase-sensitive and 
phase-insensitive correlations, and is valid in the quasi- 
monochromatic regime as well as in the broadband 
regime. Furthermore, to observe this effect with clas- 
sical fields, it is preferable to utilize narrowband sources. 
On the other hand, the marginal improvement in resolu- 
tion observed in Sec. lIVl is a strictly broadband effect that 
manifests itself in complex-stationary phase-sensitive and 
phase-insensitive correlation functions. The two experi- 
ments capitalize on different properties of phase-sensitive 
coherence and therefore they are not experiments demon- 
strating equivalent physical principles. 

In conclusion, we have presented a unified Gaussian- 
state analysis of two transverse imaging configurations, 
one that images the far-field diffraction pattern of a 
source-plane transmission mask, and one that performs 
thin-lens imaging of that source-plane transmission mask. 
We have shown that the far-field diffraction patterns 
obtained with classical phase-sensitive Gaussian-state 
light and nonclassical Gaussian-state light with low- 
brightness — such as the biphoton — differ only in con- 
trast, viz., the fringe compression is a classical phe- 
nomenon owing to the far-field diffraction of phase- 
sensitive coherence. In the second experiment, we 
have demonstrated that the cross-frequency coupling 
in complex-stationary broadband phase-sensitive light — 
whether classical or quantum — leads to a slightly nar- 
rower point-spread function than that obtained with 
quasimonochromatic phase-sensitive or phase-insensitive 
light. However, because of the enormous bandwidth that 
is necessary to observe any appreciable change in the 
point-spread function, contrary to what is stated in pT| . 
there is no practical advantage to be gained from broad- 
band operation in this image acquisition configuration. 
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